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We report on a proposal to change the group velocity of a small Bose Einstein Condensate (BEC) 
upon collision with another BEC in analogy to slowing of light passing through dispersive media. 
We make use of ultracold collisions near a magnetic Feshbach resonance, which gives rise to a 
sharp variation in scattering length with collision energy and thereby changes the group velocity. A 
generalized Gross-Pitaveskii equation is derived for a small BEC moving through a larger stationary 
BEC. We denote the two condensates by laser and medium BEC, respectively, to highlight the 
analogy to a laser pulse travelling through a medium. We derive an expression for the group 
velocity in a homogeneous medium as well as for the difference in distance, 8, covered by the laser 
BEC in the presence and absence of a finite-sized medium BEC with a Thomas-Fermi density 
distribution. For a medium and laser of the same isotopic species, the shift S has an upper bound 
of twice the Thomas-Fermi radius of the medium. For typical narrow Feshbach resonances and a 
medium with number density 10 14 cm -3 up to 85% of the upper bound can be achieved, making 
the effect experimentally observable. We also derive constraints on the experimental realization of 
our proposal. 



Over the last two decades significant advances have 
been made to replicate linear and non-linear optical phe- 
nomena with matter waves, creating the field of matter- 
wave optics. For example, atom lasers [TJ[2] are sources of 
coherent ultra-cold atoms generated by extracting atoms 
from a Bose-Einstein condensate (BEC). The coherence 
of BECs was demonstrated by interfering two conden- 
sates [3]. Atomic mirrors and beam splitters have also 
been realized [4]. Recently, the matter- wave equivalent 
of meta-materials (media with negative refractive index) 
has been proposed [5] . The analog of nonlinear four- wave 
mixing has been demonstrated using atom lasers [6j[7]. In 
these experiments three BECs with phase-matched rela- 
tive momenta generated a fourth beam. 

We present a proposal to slow a BEC while propa- 
gating through another BEC near a magnetic Feshbach 
resonance in analogy to slowing of light in dispersive me- 
dia. Slowing of light occurs when the refractive index of 
a medium varies sharply with photon frequency. Using 
electromagnetically induced transparency [8], slow light 
has been observed with a Bose-Einstein condensate [9|fTQ] 
and a hot Rb gas [11] acting as the medium. 

Magnetic Feshbach resonances are a tool with which to 
manipulate the interaction between ultracold atoms [12] . 
They are used for creating ultracold molecules, molecu- 
lar condensates , and in the BEC-BCS crossover in Fermi 
gases p~3j [14] . Feshbach resonances play an essential role 
in condensing 133 Cs, 85 Rb and 39 K [T5UT7] . Cooling ul- 
tracold atoms using Feshbach resonances has been pro- 
posed [18]. Collisions can also be tuned using optical 
Feshbach resonances [19-21], as their width can be dy- 
namically varied with a laser. They, however, tend to 
suffer from losses due to spontaneous emission. 

Figure [I] shows our proposal of slowing a condensate 
travelling through a large stationary BEC near a mag- 
netic Feshbach resonance. These two BECs play the role 



of a laser and a medium, respectively. The laser BEC 
contains atoms of mass m and has a diameter t\,. It 
propagates with wavevector k of magnitude ko ^> 1/^l- 
Its average velocity is vo = fiko / m and the kinetic energy 
per particle is Eq = h 2 ^/ (2m). (Here h is the reduced 
Planck's constant.) It is incident on a stationary BEC 
of size 4, such that t\, <C £m (in all spatial directions). 
Figures [TJd and c show two distinct cases of collisions be- 
tween condensates. In Fig.[T]3 elastic scattering out of the 
two condensates is significant. These losses have limited 
four- wave mixing experiments [22 . On the other hand, it 
allows detection of d-wave shape resonances in collisions 
between BECs [23] and thermal gases [H]. In Fig. [J 
the collision takes place at the lossless point, where the 
scattering cross section is small due to interference of the 
background and resonance scattering amplitudes. As we 
will derive the laser BEC then slows down. The slowing 
can be quantified by the difference in distance trav- 
elled by the laser BEC in the presence and absence of 
the medium. 

We first review the two-body physics of a Feshbach res- 
onance. Subsequently, we derive the equation of motion 
for a laser BEC travelling through a medium. We then 
derive the group velocity in a homogeneous medium fol- 
lowed by an estimate for 5 for an inhomogeneous medium 
whose density is described by a Thomas- Fermi profile. 
Finally, restrictions on an experimental realization are 
given. 

Collisions between ultracold atoms are dominated by 
s- wave scattering and the scattering amplitude f(E) only 
depends on the relative collision energy, E = h 2 k 2 /(2ja). 
Then around an isolated magnetic Feshbach resonance 
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FIG. 1: A schematic of a slow-atom experiment with colliding 
condensates. The left and right hand side show images of the 
condensates before and after the collision, respectively. Panel 
a) simply shows the laser BEC moving in free space with 
kinetic energy Eq and wave vector ko. Panels b) and c) show 
two cases of a laser BEC propagating through the medium 
BEC. Panel b) shows a case with large elastic scattering losses 
indicated by the halo of scattered atoms . Finally, panel c) 
depicts the collision near a Feshbach resonance where there 
is negligible scattering loss. This occurs when the scattering 
amplitude is zero. The distance delay S of the laser BEC is 
also indicated. 



where \i is the reduced mass and k is the magnitude 
of relative wave vector. The background scattering am- 
plitude fbg(E) = e^ bg sin(5b g )/fc with phase shift #b g . 
To a good approximation f^ g (E) = — «bg/(l + ikabg), 
and #b g = — kcLbg, where <2b g is the background scatter- 
ing length. The resonance width hT(E) = 2/cab g r in 
the threshold limit k 0. The energy-independent re- 
duced width Tq = /i res A, where /i res is the difference 
between the magnetic moments of the resonance state 
and the asymptotically free atoms and A is the mag- 
netic width of the resonance. The resonance energy is 
E res = ^ res (B — -Bq), where B is the magnetic field and 
Bq is the resonant field. The scattering amplitude f(E) 
satisfies the optical theorem [25] . 

Figure [2] shows f(E) near a Feshbach resonance as a 
function of collision energy E. The resonance occurs at a 
finite collision energy and f(E) approaches fb g (E) away 
from E res . The imaginary part of f(E) is related to the 
total cross section cr(E) and thus to the fraction of scat- 
tered atoms. In fact, cr(E) = Airlm f(E)/k from the 
optical theorem [25]. On resonance Im(/) is maximal 
and w 1/fc. More importantly, there exists a collision 
energy at which f(E) = due to an interference be- 
tween the background and resonance scattering ampli- 
tudes. This collision energy, which to good approxima- 
tion is E res — T(E res )/2 for positive <2b g , is the lossless 
optimal point mentioned in Fig. [I] 

We now describe the many-body physics of collid- 
ing BECs, assuming that both BECs contain the same 
atomic species. Then the collision energy E « Eq/2. 
Their dynamics is well described by the time evolution 
of the order parameter \£(x, t), the expectation value of 
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FIG. 2: (Color online) Real and imaginary part of scattering 
amplitudes, / and /b g , as a function of collision energy for 
a narrow 87 Rb Feshbach resonance at Bo = 68.54 mT and 
A = 0.6 fiT, with B - B tuned at 0.05 mT. Also shown is 
the lossless point where lmf(E) = 0. Values obtained from 
|12j . lao = 0.0529 nm, and k& is the Boltzmann constant. 



the annihilation operator ^(x,t) in the Heisenberg pic- 
ture. For a BEC at rest its evolution is well described 
by the Gross-Pitaevskii (GP) equation, derived for an 
energy-independent and real scattering amplitude. Both 
assumptions are invalid near a Feshbach resonance at fi- 
nite energy. 

Our starting point is Eq. 38 of Ref. [27 obtained using 
a cumulant expansion. It includes the time and energy 
dependence of the two-body scattering and is given by 
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where to is the initial time, Hib — —h 2 V 2 /(2m) + V(x), 
the single-particle Hamiltonian, and V(x) is the external 
potential. The operator T 2 b is the two-body T-matrix in 
the time domain and the integrals over y^ for i = 1,2,3 
are in coordinate space. 

The momenta of the medium and the laser BEC are 
non-overlapping. Hence, the wave function \P(x,t) is 
the sum of orthogonal medium and laser wavefunctions, 
^m(x, t) and t), respectively. Since the laser BEC 

is small in size and number density, we linearize Eq. ([2| 
for the laser condensate assuming phase-matching condi- 
tions. After rearranging terms, we find 

*ft^*L(x,t) = J ffiB* L ( x ,t)+2 f Y[d yi ** u (y 3 ,t) 

i 
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x / dt x (x,y 3 |T 2B (Mi) |yi,y2> 

x* L (yi,*i)*M(y2,*i), (3) 

which is nonlocal in both space and time and ignores 
the effect of the laser on the evolution of the medium. 
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Hence, the medium BEC is described by the "energy- 
independent" GP equation. 

A more insightful expression can be obtained when we 
approximate the integrands in Eq. ([3| by power series in 
derivatives evaluated at x and t. First, we realize that 
^2b(Mi) on ly depends on t — t\ and is peaked around 
t — t\ = 0. Assuming that the wave functions vary slowly 
in time, the lower limit of the integral over time can be 
extended to — oo. Next, we note for any pair of functions 
g(t) and h(t) 
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where h(z) = J dte lzt h(t) is the Fourier transform of h(t) 
and h(id/dt) = ^ n d n h/dz n \ z=0 (id/dt) n /n\. Using 
Eq. (4) with h(r) = (.| T 2 b (t) |-), the interaction term in 
Eq. (3) reduces to 

2 J Y[dy^* M (ys,t) (5) 

i 

x (x, y 3 1 T 2B {ihd/dt) I y x , y 2 > * l (y 1 , *) *m (y 2 , *) , 

where the T-matrix T 2 b(z) is now in the energy domain 
(dropping the ~ notation for simplicity) and the time 
derivatives only act on ^(yi, £)^M(y2, t). 

The T-matrix in coordinate space can be evaluated by 
transforming to the momentum representation. For s- 
wave scattering the dependence on the relative momenta 
can be neglected. That is, to a good approximation the 
T-matrix in momentum representation is [27] 
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where the J-function reflects total momentum conserva- 
tion, f(E) is the scattering amplitude of Eq. 0, and 
M = 2m. 

Inserting the momentum representation of T 2 b into 
Eq. ([5| and noting formally that 0(y,t) = exp[i(y — x) • 
zV]0(x, t), the Taylor expansion of 4>{y,t) around posi- 
tion x, the interaction term becomes 
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where Q — k4 • x + IC3 • V3 — k 2 • y 2 — ki • y 1 . Performing 
all integrations, we find a local equation of motion for 
^ L (x,£). That is, 

8 Airfr 2 
zfi-* L (x,t) = fTi B * L (x,t) - -— **,(x,t) (8) 

x/ (ihd/dt + h 2 V 2 /2M) * L (x,t)* M (x,t) • 
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FIG. 3: (Color online) Distance delay S of the laser BEC 
normalized by the Thomas- Fermi radius £m of the medium 
BEC as a function of dimensionless parameter /3 = Ub g /To, 
where Uh g and To are defined in the text. The delay for 
selected resonances assuming a peak number density of the 
medium of um = 10 14 cm -3 is shown by colored markers. 
The inset shows the group velocity v g of the laser BEC in a 
homogeneous medium BEC as a function of /3. Here vo is the 
free space velocity of the laser BEC. Markers indicate v g for 
the same selected resonances and um as in the main figure. 



Since t\, <C £m and the spread in the collision energy 
is much smaller than T(Eq/2), it is sufficient to expand 
f(z) to first order around z = Eo/2, the average relative 
collision energy, and derivatives with respect to ( x , t) 
can be neglected. For our homonuclear system the time 
evolution of the laser condensate is then given by 



-V 2 + F mf (x) 



2m* (x) 
+ F deriv (x) * L (x,t), (9) 



where m*(x) = m[l + 2a(x)]/[l + a(x)] is the 
position-dependent effective mass and a(x) = 
(47rh 2 /m)\^ M (^)\ 2 df(z)/dz with df(z)/dz evalu- 
ated at z = Eo/2. The "mean-field" potential 
VW(x) = [y(x)-(87r?i 2 /m)|vl/ M (x)| 2 /(i? /2)]/[l+2a(x)] 
contains the external potential and a potential propor- 
tional to the scattering amplitude and medium density. 
The latter contribution is analogous to the interaction 
potential in the GP equation except that the scattering 
amplitude is evaluated at non-zero energy Eo/2. Finally, 
the potential Vderiv(x) = Eoa(x)/(l + 2a(x)). The 
factor 1 + 2a (x), appearing throughout, results from the 
ihd/dt argument of the scattering amplitude. 

The operator acting on ^l(x) on the right-hand side 
of Eq. ([9| is not Hermitian as the scattering amplitude 
is complex valued. In fact, the non-Hermicity leads to 
atom loss out of both condensates, shown in Fig. (TJd as 
the halo. For a medium number density wm, the loss rate 
out of the laser condensate is umVqcf. Consequently, at 
resonance, where a ~ 87r/(&o/2) 2 , the fraction of atoms 
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remaining in the laser condensate after the collision is 
« exp (— 87rnM^M/(^o/2) 2 ). For typical values of um and 
£m almost all the laser atoms are lost at resonance. 

For our proposal we need to minimize these losses. We 
can use the lossless point where f(z) = 0, indicated in 
Fig. [2| and the total cross section is zero. The effective 
mass and the potentials in Eq. ([9| are then real. In fact, 
m*(x) > m and df/dz = abg/IV 

At the lossless point the simplest case to analyse is 
that of a homogeneous medium and V(x) = 0. The po- 
tential V m { (x) vanishes and the effective mass is uniform. 
Transforming Eq. (|9| to momentum space, we find that 
the propagation or group velocity of the laser BEC is 



v g (ko) 



loss less 771* 



v 



1 + /3 
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where the dimensionless quantity f3 = U^/Tq > and 
Uhg = (Snh 2 /m)a Dg nM is the background mean-field in- 
teraction energy. The inset of Fig. [3] shows the group 
velocity as a function of j3. The group velocity lies be- 
tween v /2 and v and approaches v /2 when f3 — ^ oo. 

We now turn to propagation through an inhomoge- 
neous medium, but still with V(x) = 0. The assumption 
£l <C £m implies that the density variation of the medium 
orthogonal to the laser propagation direction is negligible 
and we only need to treat propagation along ko passing 
through the center of the medium BEC. 

For simplicity the density profile of the untrapped 
medium is given by |^m(x)| 2 = ^m(1 — x 2 /^m)> using 
the Thomas-Fermi approximation and assume that ex- 
pansion of the medium can be neglected. Here um is the 
peak number density and £m is the Thomas- Fermi radius 
of the medium. We assume that Vd er iv( x ) <^ Eq/2 for all 
x and, hence, can apply the Wentzel-Kramers-Brillouin 
(WKB) approximation to estimate S. We find 



arctanh 
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and the dimensionless quantity f3 = U^/Yq is evaluated 
at the peak number density tim- Figure [3] shows S as a 
function of (3. The maximum 5 that can be attained by 
the laser is 2£m for f3 — » oo. 

There are several constraints on the realization of 
the proposal. Firstly, we have £\, <C £m- Secondly, 
scattering is s-wave dominated, so that /co&bg « 1 or 
E /2 <C h 2 /(2mal g ) = E hg , the Wi gner threshold limit. 
Thirdly, by solving for lmf(z) = for positive a Dg , we 
find that the requirement r < E res must hold. Fourthly, 
the energy window around the lossless point, where Im/ 
is small, is on the order of T(E res ). Consequently, the 
energy width of the laser BEC, AE L = H 2 /(2m£l) , 
must satisfy AEl <C T(E res ). In other words £\, ^> 



ft 2 /(mab g r ) = £™ in . Finally, we require resonances for 
which 5 is comparable or larger than the size of the laser 
BEC. Since £\, <C £m-> we have j3 = U^/Yq is at least of 
order one. 

TABLE I: Resonance parameters, experimental constraints, 
and spatial delay for nine Feshbach resonances. The first five 
columns specify the Feshbach resonance. The columns are 
the atomic species, magnetic resonance position .Bo, magnetic 
width A, reduced width To, and Wigner-threshold limit i^bg- 
The sixth column gives the minimum size, i™ 111 , of the laser 
BEC. The last column is the shift 5 in units of the radius of the 
medium £m, assuming a peak medium density of um = 10 14 
cm -3 . Parameters obtained from 1121. 



Atom 


Bo 


A 


T /k B 


E hg /k B 




S/£u 




(mT) 


(mT) 


(jjlK) 


(jjlK) 


(fim) 




23 Na 


119.5 


-0.14 


14 


1900 


0.45 


0.15 


a 


90.7 


0.10 


260 


1900 


0.025 


0.0091 


a 


85.3 


2.5 x 10" 4 


0.64 


1900 


9.8 


1.20 


87 Rb 


100.74 


0.021 


39 


200 


0.027 


0.025 


II 


91.17 


1.3 x 10" 4 


0.24 


200 


4.4 


1.25 


II 


68.54 


6 x 10" 4 


0.54 


200 


1.9 


0.89 


II 


40.62 


4 x 10" 5 


0.054 


200 


19 


1.7 


II 


9.13 


1.5 x 10" 3 


2.0 


200 


0.52 


0.38 


52 Cr 


49.99 


0.008 


22 


290 


0.076 


0.078 



Table [T] gives a non-exhaustive list of narrow reso- 
nances, which satisfy the constraints. For four of these 
resonances the expected S is shown in Fig. [3] assum- 
ing a peak density of um = 10 14 cm -3 . If we assume 
£l/£m ~ 0-1 then S ranges from O.I^l to 20^l fc> r the res- 
onances in Table [l[ For the selected density the chromium 
resonance is only a marginal candidate for slowing exper- 
iments. However, larger um are achievable. 

In conclusion we have shown that collisions in the 
presence of a magnetic Feshbach resonance can lead 
to slowing of a laser BEC as it propagates through a 
large medium BEC. The slowing is a consequence of 
the collision-energy dependence of the scattering ampli- 
tude near the resonance. Based on a generalized Gross- 
Pitaevskii equation, we predict a maximal reduction of 
the group velocity by a factor of two and suggest that the 
experiment be performed at a magnetic field where the 
elastic scattering is zero. Such a field always exists near 
a magnetic Feshbach resonance. For finite-sized conden- 
sates slowing can be observed by measuring the spatial 
delay of the laser BEC, which can not exceed twice the 
Thomas-Fermi radius of the medium. We show that for 
narrow resonances this signal is expected to be measur- 
able. 
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